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Abstract 

We consider the highly anisotropic ferromagnetic spin 1/2 Heisenberg chain with pe- 
riodic boundary conditions. In each sector of constant total z component of the spin, 
we develop convergent expansions for the lowest band of eigenvalues and eigenfunctions. 
These eigenstates describe droplet states in which the spins essentially form a single linear 
droplet which can move. Our results also give a convergent expansion for the dispersion 
relation, i.e., the energy of the droplet as a function of its momentum. The methods used 
are from [4] and [5], and this short paper should serve as a pedagogic introduction to 
those papers. 
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1 Introduction 



We consider the spin 1/2 ferromagnetic XXZ Heisenberg chain. With periodic boundary con- 
ditions its Hamiltonian is 

H = - XK^i + + (!) 

We will always assume that e is small, so we are in the Ising-like regime. This Hamiltonian 
commutes with the total z component of the spin, so one can study it in a sector with a fixed 
total z component of spin. The lowest energy states in each such sector are droplet states in 
which the down spins essentially form a single droplet. This droplet can move around, so the 
lowest part of the spectrum in each sector should be a band of continuous spectrum. Much 
of this picture was proved by Nachtergaele and Starr [14]. In particular, they showed that 
the eigenvalues of these droplet states lie in an interval whose width is of order e m where m 
is the number of spins in the droplet, and there is a gap between this band and the rest of 
the spectrum. They gave an explicit description of the subspace spanned by these droplet 
eigenstates which becomes exact in the limit that the size of the droplet goes to infinity. 

We will show that methods from [4] and [5] can be used to construct convergent expansions 
for the eigenstates in this lowest band in each sector with a fixed number of down spins. Because 
the translation-invariant ground states of the ferromagnetic XXZ chain are trivial, the estimates 
in this paper are considerably simpler than those in [4] and [5]. Thus this short paper should 
serve as a pedagogic introduction to those papers. At first glance it may be surprising that 
one can construct convergent expansions for eigenstates that will become part of continuous 
spectrum in the limit that the length of the chain goes to infinity. We avoid this problem 
by using the fact that with periodic boundary conditions the eigenfunctions have a definite 
momentum. Within the subspace of a given momentum the eigenvalue we want to study is 
isolated, and a convergent expansion is possible. The results of [14] are valid for |e| < 1, while 
our results are valid only for small e. So our results have a smaller range of applicability, but 
give more detail about the droplet eigenstates and eigenvalues. 

Before we consider the droplet states, we will first consider the chain with the periodic 
boundary conditions replaced by boundary conditions that favor the spin to be up at one end 
of the chain and down at the other end of the chain. Such boundary conditions force a kink or 
interface into the chain. For a generic Hamiltonian one would not expect this interface to be 
stable. Roughly speaking it would be equally likely to be anywhere along the chain. However, 
for the model considered here the conservation of total z component of spin essentially fixes 
the location of the kink. The stability of this interface was proved independently by Alcaraz, 
Salinas and Wreszinski [1] and Gottstein and Werner [6]. Alternative proofs of the stability 
of this interface were given in [3], by using the path integral representation of interface states, 
and in [2]. Other work on these interface states includes [7], [8], [9], [10], [11], and [12]. Our 
discussion of these interface states does not contain any new results and is intended as an 
introduction to the methods used for developing the expansion of the droplet states. 
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We will always work on finite chains with estimates independent of the length of the chain. 
Of course, one obtains continuous spectra only in the infinite length limit. It would be inter- 
esting to show that in the infinite length limit, defined for example via the GNS construction, 
there is indeed a band of continuous spectrum corresponding to the motion of the droplets. 

This paper only concerns the one- dimensional XXZ ferro magnet. There are many interesting 
results on interfaces in the higher dimensional model (see [13] for a review), but we are not 
aware of any results on droplet states in the higher dimensional model. In two dimensions, if 
m, the number of spins in the droplet, is a perfect square, then the ground states of the e = 
Hamiltonian are just droplet states in which the droplet is a square. So the e = ground state 
is unique up to translations. In these sectors it may be possible to use the methods of this 
paper to contruct the droplet eigenstates. (A similar remark applies in higher dimensions.) For 
other sectors the e = ground states are typically more degenerate, and the methods of this 
paper would need to be combined with some form of degenerate perturbation theory. 

2 Kink states 

In this section we construct convergent expansions for the "kink states" which contain a single 
domain wall. These states have been studied extensively, and there are no new results in this 
section. This section serves as an introduction to our expansions for the droplet states in the 
next section. The expansion for the kink states is similar to that for the droplet states, but 
simpler since one does not have to deal with the freedom of the droplets to move about. 

We consider the ferromagnetic XXZ Heisenberg chain on N sites with the following boundary 
conditions. 

N-l 

H= £[1- o]o) +1 - ea*a* +1 - eo)a)^\ + A(l + a{) + B(l - a z N ) (2) 

The constants A and B will be positive, so the boundary conditions favor the spin down state 
at site 1 and the spin up state at site N. If we consider the sector with m down spins, then the 
above Hamiltonian with e = has a unique ground state given by putting the to down spins 
at sites Horn and the up spins at the other sites. For a set X of sites we define \X) to be the 
state with down spins at the sites in X and up spins everywhere else. So the e = ground state 
is |{1, 2, • • • , to}). The set {1, 2, • • • , to} will appear throughout our equations, so we denote it 
by 

M = {1,2,---,to}. (3) 

When e is not zero, the ground state in the sector with to down spins will be a perturbation of 
\M). The coefficient of \Y) should get smaller as Y gets "farther" from M. Thus it is natural 
to write Y as a modification of M. For sets of sites X and Y we denote by XAY the set of 
sites which belong to exactly one of X and Y. So 

XAY = (X U Y) \ (X n Y) (4) 
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Any state in the sector with to down spins can be written as 

= ]T e(A) | A" AM) (5) 

for some coefficients e(X). The constraint X : to means that we only sum over sets X such 
that | A AM) has to down spins. This means that AAM has to elements, or equivalently that 
X (1 M and X n M c have the same number of elements. For example, X could be {to, to + 1} 
in which case |MAA) is the state that occurs at first order in perturbation theory. 

When e = 0, \XAM) is an eigenstate for all X. We let X(X) denote its eigenvalue. It 
equals 2 times the number of nearest neighbor pairs of sites for which the spins in \XAM) are 
anti-parallel plus boundary terms. The boundary terms are +2 A if 1 e X and +2B if N e X. 
For the groundstate, X = 0, we have X(X) = 2. If we have A > 1 and 5 > 1, then for all 
other X with X : to we have A(A) > 6. (Note that there are other X with A(A) = 2, e.g., 
X = {to+ 1}, but they belong to sectors with a different number of down spins.) If A and B are 
positive, but possibly less than 1, then the lower bound of 6 must be replaced by the minimum 
of 6 and 2 + 2 A + 2B. The following proof applies in this case with slight modifications. To 
keep things as simple as possible we will assume A and B are both at least 1. (In the more 
general case, just how small e would need to be would depend on A and B.) 

We now consider the Hamiltonian with nonzero e. The term cr? cr? +1 + a|o"| +1 can be written 
as (j|(jj +1 (l — <t|<7| + i). So if the spins at j and j + 1 are parallel, it gives zero and if they are 
anti-parallel it flips these two spins and multiplies the result by 2. We let d(XAM) denote the 
set of bonds < j, j + 1 > such that the spins at j and j + 1 in \XAM) are different. One can 
think of d(XAM) as the set of domain walls in the spin configuration. Then 

H* = £ X(X)e(X) \XAM) + 2e £ e(X) £ j + 1}AXAM) (6) 

X:m X:m <j,jf+l>ed(XAM) 

Note that < j, j + 1 >G <9(XAM) if and only if < j,j + 1 >G d({j,j + 1}AXAM). So we can 
do a change of variables X — > XA{j, j + 1} in the last term and rewrite this as 

H* = £ ApQepQ |XAM) + 2 e £ £ e(XA{j, j + 1}) |XAM) (7) 

X:m X:m <j j+l>e9(XAM) 

When e = 0, the ground state energy is 2, so we will write the ground state energy for nonzero 
e in the form 2 + E, i.e, we look for a solution of Hty = (2 + E)ty. Equating the coefficients of 
| A AM), we have for A with A : to, 

A(A)e(A) + 2e ]T e(AA{j, j + 1}) = (2 + £)e(A) (8) 

<j,j+l>€d(XAM) 

We take e(0) = 1. (At first we cannot be sure this is possible since e(0) could be zero for 
the groundstate. But if we can succeed in constructing an eigenstate with this condition and 
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show it is the ground state, then this will justify the assumption.) The unknowns in the above 
equation are e(X) for non-empty X with X : m and E. For X = the above equation gives 

E = 2ee({m,m + 1}) (9) 

and for non-empty X with X : m the above equation can be written as 

1 



e(X) 



X(X) -2 



-2e E e(XA{j,j + l}) + £e(X) 
<i,i+i>e9(XAM) 



(10) 



We think of equations (9) and (10) as a fixed point equation for the unknowns E and e(X), 
where X ranges over non-empty sets with X : m. We use e to denote this collection of unknowns 
and write equations (9) and (10) together as F(e) = e. We define a norm on the space of e's 
by 

INI = \E\+ E (A(X)-2)|e(X)| (11) 

X:m,X^V) 

We will show that F is a contraction in a small ball about the origin and maps this ball back 
into itself. The existence of a solution then follows by the contraction mapping theorem. 

At first glance it seems that F(0) = 0, but it is not. This is because e(0) can appear in the 
right side of (10). This happens when X = {m,m + 1} and j = m. This produces a "constant" 
term in F(e), i.e., a term which does not depend on e. Denoting the terms in F(0) by cq{X) 
we see that e ({m,m + 1}) = — e/2, and the other e (X) and E are all zero. We note, for later 
use, that ||F(0)|| = 2e. 

To show that F is a contraction we have 

||F(e) -F(e')|| < 2e|e({m, m + 1}) - e'({m, m + 1})| 

+ 2e E E |e(XA{j,j + l})-e'(XA{j,j + l})| 

+ E |£e(X)-£V(X)| (12) 

We split this bound into two parts. The first term in the right side of the above is just what 
the second term would give when X = 0. So the first two terms of the three terms in the right 
side of the above are 

= 26 E E \e(XA{j,j + l})-e'(XA{j,j + l})\ 

X:m <j,j+l>£d(XAM) 

= 2eJ2\(X)\e(X)-e'(X)\ (13) 

X:m 

where the last equality follows by a change of variables, X — > XA{j,j + 1}. For X = 0, 
\e(X)-e'(X)\ = 0. For all other X with X : m we have \{X) > 6 and so X(X) < |(A(X) -2). 
Thus the above is 

<3e||e-e'|| (14) 
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The second part of (12) is 

E \Ee(X)-E'e'(X)\ < £ [ |£| |epT) - e'pT)| + \E - E'\ \e\X)\ ] 

X:m,X^ X:m,X^9 

< \ E (X(X)-2)[\E\\e(X)-e'(X)\ + \E-E'\\e'(X)\] 

< J||e-e'|| max{||e||, ||e'||} (15) 

where we have used X(X) — 2 > 4. 

Putting together these two bounds, we have shown for e,e' with norm less than 8 that 

||F(e)-F(e')||<(36 + ^max{||e||,||e'||})||e-e'||<(3e+^||e-e'|| (16) 

Thus F is a contraction on the ball of radius 5 about the origin if 5 and e are small enough. 
To see that it maps the ball of radius 5 about the origin back into itself, we use 

1 1^)1 1 < + \\F(e)-F(0)\\ < 26+ (36+ 6 -)\\e - 0|| < 5 (17) 

if ||e|| < 5 and e is sufficiently small. 

By using a stronger norm we can derive decay properties for the coefficients e(X). Let K 
be a large positive constant. Define 

||e|| = |£|(K| e |)- 2 + E (A(X)-2)|e(X)|(K| e |)-W (18) 

X:m,X^9 

We will assume that e is small enough that K\e\ < 1. w(X) is a positive integer which is 
the lowest order in perturbation theory at which e(X) gets a nonzero contribution. More 
precisely, we consider all sequences of sets X , Xi, ■ ■ ■ , X n such that X = 0, X n = X and for 
% — 1, 2, • • • , n, Xi — Xj_iA{j, j + 1} for some j G 9(Xj_iAM). Then w(X) is the smallest n 
for which such a sequence exists. 

The previous estimates that proved the existence of a fixed point in our original norm can 
be repeated with this new norm. One finds that the estimates continue to hold provided w(X) 
satisfies the following three properties. 

w({m, m + 1}) = 1 
w(XA{j,j + l})<w(X) + l 
w(X) - 2 < w(X) 

(19) 

The first two properties follow easily from the definition of w(X), while the third is trivial. 
Note that the existence of a solution to the fixed point equation in this stronger norm implies 
that the coefficients e(X) decay at least as fast as (K\e\) w<yX \ 
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The expansion that we have developed can now be used to study a variety of properties of 
the kink states. For example, the localization of the kink near the site m follows from these 
estimates. One can consider the dependence of the ground state energy on the sector (the 
choice of m). Most studies of the kink states used particular values of A and B for which this 
energy is independent of the sector. Bach and Macris [2] considered more general boundary 
conditions and showed that the difference between the ground state energies in different sectors 
were exponentially small in the length of the chain provided the kink is not near the boundaries. 
With more work it is probably possible to rederive this result with our expansion. Our interest 
here is primarily in the droplet states, so we do not pursue this approach to the kink states any 
further. 



3 Droplet states 

Now we consider the ferromagnetic XXZ Heisenberg chain with periodic boundary conditions. 
For N sites its Hamiltonian may be taken to be 

N 

H = Et 1 - ~ " (2°) 

i=i 

Indices will always be taken to be periodic, e.g., cr z N+1 means o\. 

We consider the sector with m down spins. We continue to use the abbreviation 

M = {l,2,-..,m} (21) 

For e = the ground states in this sector are |M) and its translates. When e ^ 0, the interface 
between sites m and m + 1 and between sites 1 and N will spread out somewhat. States of the 
form | MAX) where X is small and localized near 1 and m will make up the dominant part 
of the eigenstate. We continue to denote the constraint that \XAM) is in the sector with m 
down spins by X : m. We now look for eigenstates with momentum k in the form 

N 

= £ e M E < X ) I (XA.M) + (22) 

1=1 X:m 

For a set of sites Y, we use Y + 1 to denote the translate of the set by / sites to the right. So 
Y + I = {i + I : i £ Y}. This state has momentum k in the sense that if T is the operator of 
translation by one lattice site to the right, then = e~ lk ^k- 

As before, dX is the set of bonds such that one endpoint is in X and the other endpoint 
is not in X. We define n(X) to be the number of bonds in <9(XAM). Then 2n(X) is the 
eigenvalue of \XAM) when e = 0. We have 

N 

H^ k = 2 ]T e M ]T e(X)n(X) | (X AM) + I) 

1=1 X:m 
N __ 

+2e£e^£e(X) £ \{j,j + l}A[(XAM) + l}) (23) 

1=1 X:m <j,j+l>ed({XAM)+l) 
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Using the change of variables j — > j + / in the second term, this is 

N 

= 2 E e m E e(X)n(X) | (X AM) + I) 

1=1 X:m 

+2e E e m E E e(X) \ ({j, j + 1} AX AM) + I) (24) 

1=1 X:m <j,j+l>£d(XAM) 

In the second sum on X we do a change of variables: X — > XA{j, j + 1}. Since < j,j + 1 >G 
<9(XAM) if and only if < j,j + 1 >G d({j,j + 1}AXAM), we obtain 

N 

HV k = 2 E e ikl E e{X)n(X) \ (XAM) + I) 

1=1 X:m 
N 

+2eJ2e M E E e(XA{ J)J + l})|(XAM) + (25) 

J=l X:m <j,j+l>ed(XAM) 

The above should be equal to E(k)^ k . The eigenvalue now depends on k. We write E(k) 
as a Fourier series. Since it equals 2 when e = 0, we take the series in the form 

N 

E(k) = 2 + E e s e ihs (26) 

s=l 

So 

V N 

(E(k)-2)* k = e e * e * fes E eiM E e ( x ) I {XAM) + o 

s=l J=l X:m 

V V 

= E E e « eiW E e ( X ) I (XAM) +l-s) (27) 

s=l /=1 X:m 

where we have used the change of variables I — > I — s. Define Y by 

(XAM) + 1 - s = (YAM) + / (28) 
Then solving for X in terms of Y we find 

X = (Y + s)A(M + s)AM (29) 

So 

N N 

(E(k) - 2)^ k = E E e *e ikl E <( Y + s ) A ( M + s ) AM ) I ( rAM ) + (3°) 



s=l 1=1 X:m 
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We now multiply both of (25) and (30) by e %kn and sum on k. This yields 

2 e ( x )( n ( x ) - 2 ) \(XAM) + n) 

X:m 

+2e J2 E e(XA{j,j + 1}) \{XAM)+n) 

X:m jed(XAM) 

N 

= E e (( Y + s)A(M + s)AM) \(YAM) + n) (31) 

s=l X:m 

Thus for sets X such that X : m we have 

2(n(X)-2)e(X) + 2e £ e(XA{j,j + l}) 

jed(XAM) 

N 

= Y / e s e((X + s)A(M + s)AM) (32) 

The smallest n(X) can be is 2. It attains this value for the empty set and for X of the form 
MA(M + n) for some n. (These are the sets for which XAM is just a translate of M.) We take 
e(0) = 1 and e(X) = for all other X with n(X) = 2. It is not clear at first that we can do this, 
but if we can succeed in constructing the eigenf unctions under this condition that will show 
that these conditions can be imposed. We note that since we now have N eigenfunctions, there 
are N degrees of freedom corresponding to their normalizations. The conditions we impose can 
be thought of as fixing these normalizations. 

If X = MA(M — n), then 

N 

Y,e s e((X + s)A(M + s)AM)=e n (33) 

s=l 

So for these X, eq. (32) becomes 

e n = 2e ]T e(MA(M-fi)AO,i + l}) (34) 

jed(M-n) 

For X with n(X) > 2, we solve eq. (32) for e(X) : 



jed(XAM) 



+ 2 n { x)-2 ? x 6s e((X + S)A(M + S) AM) (35) 

We use e to denote the collection of variables e(X) for X such that X : m and n(X) > 2 
and the variables e„. These are the unknowns in (34) and (35). The right sides of (34) and 
(35) define a function F(e), and together these two equations can be written as the fixed point 
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equation F(e) = e. It is important to ask if the argument of any of the e( ) in the right sides 
of these equations can be the empty set or a set of the form MA(M + s). It is easy to check 
this does not happen in (34), except for the trivial case of m — 1. And it does not happen in 
the second term in (35) thanks to the constraint n(X) > 2. But it can happen in the first term 
in (35). Since e(MA(M + s)) = for s ^ 0, these terms drop out of (35). And since e(0) = 1, 
when X = {m, m+l}orX = {1, N} we get a contribution of — e/2. In particular, F(0) is not 
the zero vector in our Banach space. 

We first show the fixed point equation has a solution using the norm 

N __ 

||e|| = £|e n |+2]r|e(X)|(n(X)-2) (36) 

n=l X:m 

We will show that with the above norm, F is a contraction in a small ball about the origin. 
We have 

\\F(e) - F(e')\\ < 2e]T ]T \e(MA(M - n)A{j,j + 1}) - e'(MA(M - n)A{j,j + 1})| 

n jed(M-n) 

+ 2e J2 E \e(XA{jJ + l})-e'(XA{jJ + l})\ 

X:n(X)>2,X:m jed(XAM) 

+ E El e - e (( X + s ) A ( M + s ) AM ) -e' s e\(X + s)A(M + s)AM)\ (37) 

X:n(X)>2,X:m s 

The terms in the first sum in the above are the terms one would get from the second sum for 
X with n(X) = 2. So together these two sums are 

= 26 £ ]T \e(XA{j,j + l})-e'(XA{j,j + l})\ 

X:m jed(XAM) 

+ 2eJ2 £ \e(X)-e'(X)\ = 2eJ2n(X)\e(X)-e>(X)\ 

X:m jed(XAM) X:m 

(38) 

where we have used the change of variables X — > XA{j, j + 1}. If n(X) = 2, then e(X) = e'(X). 
For the other X, n(X) > 4 and so n(X) < 2(n(X) — 2). So the above is < 4e||e — e'||. 

In the third sum in (37), we drop the constraint n(X) > 2 and do a change of variables 
Y = (X + s)A(M + s)AM. Then it is 

< EEW)-^'00l 

Y:m s 

< EE \e.\ HY) - e'OOl + E E K - 41 |e'(r)| 

< max{||e||,||e / ||}||e-e / || < «J||e-e'|| (39) 

Thus F is a contraction on a ball of radius 5 about the origin if 8 and e are small enough. Since 
||F(0)|| is of order e, the above estimate also shows that F maps the ball {e : ||e|| < 6} back 
into itself if 8 and e are small enough. 
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As in the previous section, we can introduce a stronger norm. Let 

l|e|| = E \e n \(K\e\)-^ + 2 £ \e(X)\(n(X) - 2){K\e\)^ (40) 

n=l X:m 

where w(X) is defined as in the previous section and w n = w(MA(M + n)). We emphasize 
that this is a natural norm in the sense that the powers w n and w(X) are the lowest order in 
perturbation theory at which the corresponding terms get nonzero contributions. We leave it to 
the reader to check that the preceding estimates go through if K\e\ < 1 and K is large enough. 
It is easy to see that w 1 — m and the other w n are even larger. Thus this stronger norm shows 
that the coefficients in the Fourier series of the dispersion relation are at least order e m . 
We summarize our results on the droplet eigenstates in a theorem. 

Theorem 1 For the chain with periodic boundary conditions and N sites, we consider the 
sector with m down spins (0 < m < N). There is a constant K > (independent of N) 
such that if K\e\ < 1, then the fixed point equation has a solution in the above norm. The N 
eigenstates and eigenvalues determined by this solution through (22) and (26) have the lowest 
eigenvalues in this sector. (Except for k = in which case it is the second lowest.) Letting 
Ex{k) denote the eigenvalue for a chain with N sites corresponding to momentum k, there are 
constants d s for s G Z such that 

oo 

lim E N {k) = 2 + V d s e iks (41) 

N^oo — ' 
s=— oo 

The Fourier coefficients d s are absolutely summable, and for s ^ the coefficients are of order 
at most e m . 

The existence of a solution to the fixed point equation implies that we have constructed iV 
eigenvalues and eigenf unctions. To see that for each k the eigenvalue is the lowest eigenvalue in 
the sector of momentum k (or second lowest if k = 0) we argue as follows. We know the claim 
is true when e = 0. For a finite chain the eigenvalues are continuous in e, and so our eigenvalue 
can cease being the lowest as e is increased only by crossing another eigenvalue. At the value of 
e where such a crossing would occur, our eigenvalue would be degenerate. One can then show 
that this implies the solution of the fixed point equation is not locally unique, contradicting 
the contraction mapping theorem. (More details may be found in [4].) 

For a finite chain the eigenvalues E^(k) are only defined by a finite set of values which 
depends on N. But we can use use (26) to extend the definition to all k and so make sense of 
the limit in (41). The existence of the limit as — > oo of the Fourier coefficients e s follows by 
standard arguments. (See [4] for similar arguments.) 
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